Weak continuous monitoring of a flux qubit using coplanar waveguide resonator 
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We study a flux qubit in a coplanar waveguide resonator by measuring transmission through the system. In 
our system with the flux qubit decoupled galvanically from the resonator, the intermediate coupling regime 
is achieved. In this regime dispersive readout is possible with weak backaction on the qubit. The detailed 
theoretical analysis and simulations give a good agreement with the experimental data and allow to make the 
qubit characterization. 



Implementation of a scalable solid state quantum system 
for either quantum information processing devices or quan- 
tum limited detectors requires high-fidelity read-out methods 
for its quantum states. For superconducting qubits the switch- 
ing of a dc SQUID or a probe junction to the normal state was 
used for the determination of the qubit states Hl-ll] . However, 
the switching to the normal state leads to strong backaction 
from the measurement device on the qubit. In order to min- 
imize detector-qubit backaction the dispersive type of mea- 
surements has been proposed Hi and implemented fll?!]- In 
the frame of this method a qubit is coupled to a resonator with 
a resonance frequency much lower than the qubit transition 
frequency. The phase of the oscillator modes was monitored 
to perform qubit measurements. 

If the qubit 's characteristic frequency is comparable to the 
eigenfrequency of the resonator, coherent interaction between 
the quantum oscillator and the artificial atom (qubit) can be re- 
alized. In the optical domain this is known as cavity quantum 
electrodynamics (cQED) which can precisely describe the en- 
tanglement of a single two-level quantum system (e.g. spin 
1/2) and a cavity field [8]. Solid state cQED based on super- 
conducting circuits has been suggested [9] and demonstrated 
using superconducting qubits coupled to a supercon- 
ducting transmission line resonator. To observe resonator- 
qubit entanglement the coherent coupling strength between 
them should exceed the rate of incoherent processes in this 
system. Quantitatively this means that g > K\y where fig is 
the coupling energy between qubit and resonator, and K is the 
resonator loss rate and y is the qubit dephasing rate (so-called 
strong coupling regime). 

Besides studying the fundamentally interesting strong cou- 
pling regime, the same architecture can be used for the read- 
out of qubits as well [12]. In the dispersive regime, when the 



qubit transition frequency is far from the resonator frequency, 
even a quantum non-demolition (QND) measurement can be 
performed provided the photon numbers in the resonator are 
low enough. The phase change in the resonator at a particu- 
lar (probe) frequency can be monitored to perform qubit state 
readout 1.13.1 . This measurement can be done in either a weak 
continuous llOllU] or pulse regime 

In the dispersive regime, when the resonator acts as a qubit 
detector, the resonator-qubit backaction can be minimized by 
decreasing the coupling g. However, in order to measure 
resonator-qubit dynamics the condition g > K should be ful- 
filled in order to avoid the performance degradation of the 
resonator as a detector. In this work, we test experimentally 
the resonator-qubit system with intermediate coupling where 
g> Khuig^y. 

The resonator was fabricated by e-beam lithography and 
CF4 reactive-ion etching of a 200 nm thick Nb-film deposited 
on an undoped silicon substrate. The length of the resonator's 
central conductor is L = 23 mm, its width 50 jim and the gap 
between the central conductor and the ground plane 30 /im, 
which results in a wave impedance of about Z = 50 ^1 and 
a resonance frequency Odyjln ^ 2.5 GHz for the fundamen- 
tal, half-wavelength, mode. The resonator is coupled to the 
in- and output waveguides via two identical gaps of 90 /im 
in the central conductor corresponding to a coupling capaci- 
tance Ck- 6 fF. At temperatures T < 50 mK, the resonator 
quality factor Q = (Oy/ /c ^ 2 • 10^, corresponding to k/Ik ^ 
0.13 MHz determined by loading loss in the external 50 
impedance. In the middle of the resonator, the central con- 
ductor is tapered to a width of 1 /im and a length of 30 jim 
with 9 /im gap (see Fig. [Ha)), which provides the necessary 
qubit-resonator coupling. 

The Al persistent current qubit is fabricated in the central 
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FIG. 1 : (a) A micrograph of the cental part of the resonator. The cen- 
tral wire and the ground plane are tapered, (b) Electron micrograph 
of the qubit and the central line of the resonator 



(tapered) part of the resonator by using conventional two- 
angle shadow evaporation technique. It is galvanically de- 
coupled from the resonator line (Fig.[Hb)), making the qubit 
fabrication independent on the resonator and providing flexi- 
bility in design and fabrication of more complicated circuits. 
The qubit loop (size of 2.5 x 5 /im^) is inductively coupled 
to the resonator. The coupling constant g/27t = ?> MHz has 
been numerically estimated and is of the same order as the de- 
cay rates of flux qubits The loop is interrupted by 
three Josephson junctions. Two of them have a nominal size 
of 800 X 180 nm^, while the third is about 35 % smaller. 

Two cryoperm shields and one superconducting (lead) 
shield enclosed the sample in order to minimize the influence 
of external magnetic fields. The sample was thermally an- 
chored to the mixing chamber of a dilution refrigerator, pro- 
viding a base temperature of less than 10 mK. Measurements 
of the amplitude as well as the phase (p of the transmitted sig- 
nal through the qubit-resonator system as a function of applied 
power and frequency were performed by using a network an- 
alyzer. Attenuators are located at different temperature stages 
on the input line preventing thermal noise of those stages to 
reach the sample. The transmitted signal is amplified by a 
cryogenic amplifier placed at the 4.2 K stage and several room 
temperature amplifiers. An isolator (circulator with one ter- 
minated port) placed at the mixing chamber is used to protect 
the sample from backaction of the amplifier. The qubit energy 
bias is controlled via two small external coils located under- 



neath and above the sample. 

The qubit-resonator system in the natural (flux) basis is de- 
scribed by the Hamiltonian lfl5l] 



H 



(1) 



where £{^x) = 2<l>o/p(l - ^x/^o) is the bias of the flux qubit 
controlled by the external dc magnetic flux <i>x, = h/2e is 
the magnetic flux quantum and Ip is persistent current of the 
flux qubit Gx^z the Pauli matrices, A is the tunnel splitting 
of the flux qubit, and a\ a are the photon creation and anni- 
hilation operators. After transformation to the eigenbasis of 
the qubit and neglecting the small diagonal terms, the Hamil- 
tonian can be rewritten as 
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hCOa 
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(2) 



where ge = g' A/V^^ + is the normalized coupling, and 
hcOq = ve^pA^ is the energy level separation of the qubit. 

For n photons in the resonator, the qubit-resonator dynam- 
ics is completely confined to a two dimensional subspace with 
basis \g)\n) and \e)\n—l). The eigenenergies relative to the 
ground state of the Hamiltonian of Eq. dJ]) can be expressed 
in the analytical form \^ 



A^4 




(^+1), (3) 



were 5 = (0^ — ft)^ is the qubit-resonator detuning. By mea- 
suring the transmission through the resonator, the transition 
frequencies of the qubit-resonator system can be probed. If 
there is no interaction between the qubit and the resonator, the 
system should respond only at frequencies corresponding to 
the eigenenergies of the qubit hcOq and resonator hcOr. 

Fig. [2la) demonstrates the experimentally measured trans- 
mission amplitude t. The finite coupling g allows to observe 
anticrossings. The crossing points are visible at COq = (Or (dot- 
ted white lines in Fig. Ha)). The white solid lines were cal- 
culated using Eq. ^ for the lowest photon transition, which 
corresponds to the mean photon number in the resonator less 
then one. 

In contrast to strong coupling where two well-resolved 
spectral lines separated by the vacuum Rabi frequency at the 
resonance point (5 = 0) have been observed [10], we have 
detected the broadening of the transmission peak only and 
the peaks disappeared at 5 = 0. Numerical simulations show 
that this is a feature of the intermediate coupling when in- 
deed g ^y. However, in spite of the fact that the coupling 
is relatively weak, we have clearly observed the anticrossing 
features. 

To investigate the qubit-resonator dissipative dynamics in 
details, we analyze the Markovian master equation for the 
density matrix p 



(4) 
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FIG. 2: (Color online) a) Normalized transmission amplitude of the 
resonator as a function of the qubit energy bias e and the driving 
frequency The data were measured at a nominal temperature of 
the mixing chamber below 10 mK, ensured ksT ^h(0q^h(0r. The 
calculated frequencies of the lowest photon transition in the qubit- 
resonator system are depicted as white solid lines, b) Theoretical cal- 
culations. The number of photons is taken to be less than one, other 
parameters are taken from the experiment. 7 is adjusted for better 
correspondence of the transmission in vicinity of anticros sings. 



where the dissipative (Lindblad) term L = Lr -\- Lq presents 
two incoherent processes: dissipation in the resonator (photon 
decay) 



Lr = — (lapa^ — a^ap — pa^a) 



(5) 



and qubit decoherence: relaxation with rate 71 and pure de- 
phasing with rate 7^ 

Lq = ^(iG-pa^ - G^G-p -pG^G-) + ^(a,pa, - p), 

(6) 



where (7± = (Cx ± iGy)l2. 

We consider a ID (coplanar waveguide) resonator of size 
L (— L/2 <x <L/2) driven by the external voltage yin{t) = 
Vin COS (Ot from the left hand side through the pointlike cou- 
pling capacitance Cx located at x = —L/2. The excited volt- 
age field within the resonator can be presented as V(x,^) = 
-Vriae-'"^' ^a^e''^^)sm{7tx/L), where = ^fi(o/Cr and 
Cr is the resonator central line capacitance. The interaction 
Hamiltonian is Hint = CKVin{t)V{—L/2j) and the total driven 
system Hamiltonian in the rotating wave approximation in 
vicinity of the qubit-resonator resonance cOr ~ cOq and with 
small detunings 5r = CO — COr and 5q = CO — cOq is 



H > 



hSa 



G^ — hSrO^a + hgs {a^ G- 



(7) 



where h£l = CKVinVr- 

Solving the stationary master equation Eq. (p =0), 
one can calculate the observable quantities. In particular, the 
expectation value of the photon field in the resonator in the 
considered weak driving limit (with the truncated the photon 
space to two states) is 



{a) 



(8) 



where 5/ = 5^ + /k*/2, 5^ = 5q-\- iy are redefined detunings 
and 7 = 7i/2 + 70 is the total qubit dephasing. The am- 
plitude of the field at the output of the resonator (at x = 
L/2 + 0) is Vout{t) ^ -icoC^ZVr{{a)e-''^' + (a^)^^'^O- And 
using the explicit expression for the photon relaxation rate 
K ^ 2ZC^co^ /Cr, we find the transmission amplitude (defined 
as t = Vout/Vin) to be 



t = 



iK/2 



(9) 



In order to characterize the qubit energy, we carry out spec- 
troscopic measurements in the following way. Two waves are 
applied: an excitation wave, which excites the qubit transi- 
tions and thereby results in a change of the resonator transmis- 
sion, and a probing one close to COr, which enables to measure 
the transmission. Using the qubit parameters, t is calculated 
according to Eq. (|9]) and plotted in Fig. [3b) as a density plot 
versus 5^ and e. In this plot, 7 was used as an adjustable pa- 
rameter for better correspondence with the experimental data 
in Fig.ISa) and found to be 7/27r = 0.26 MHz. 

To understand the effect of the qubit state on the resonator 
transmission far away from the anticrossings (|5| gg), the 
Hamiltonian in Eq. Q can be transformed |9i] to 



H^fiico, 



ge 



COa 



5 



(10) 



The first term contains the ac Zeeman shift of the resonance 
frequency, which depends on the state of the qubit as well as 
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FIG. 3: (Color online) Transmission phase of the resonator as a 
function of the qubit excitation frequency e/h. The black line depicts 
the fitted energy level splitting of the qubit. 

on its detuning. The ac-Zeeman frequency shift is visualized 
by the color profile in Fig.[3l It is positive (negative) for 5 < 
(5 > 0). If the driving microwave field is in resonance with the 
qubit energy levels, the expectation value of becomes zero, 
resulting in disappearance of the ac Zeeman shift. Therefore, 
the green traces appear in Fig. [3] following the relation for 
energy levels separation of the qubit C0q{£). By fitting this 
curve (black solid line) the persistent current Ip = 180 nA as 
well as the gap A/h= 1.8 GHz of the qubit were determined. 
A similar effect was observed by Schuster et al with a charge 
qubit Q. 

In conclusion, we demonstrate that in the galvanically de- 
coupled resonator-flux qubit system, an intermediate cou- 
pling, i.e. g ^ 7, is possible. With the intermediate coupling, 
the main quantum mechanical features of this system are ob- 
servable: the <2C-Zeeman shift, the anticrossing between the 
single photon in the resonator and the qubit, as well as the 
level splitting of the qubit versus external bias. 

We gratefully acknowledge the financial support of the EU 
through the EuroSQIP and Solid project and the financial sup- 
port by the DFG IL150/b-l. EI acknowledges the financial 
support from the Federal Agency on Science and Innovations 
of the Russian Federation under contract N 02.740.11.5067. 



M.G. was partially supported by the Slovak Scientific Grant 
Agency Grant No. 1/0096/08, the Slovak Research and De- 
velopment Agency under the contract No. APVV-0432-07 
and No. VVCE-0058-07, ERDF OP R&D, Project CE QUTE 
ITMS 262401022, and via CE SAS QUTE. 



[1] Y. Nakamura, Y. A. Pashkin, and J. S. Tsai, Nature 398, 786 

(1999) . 

[2] C. van der Wal, A. ter Haar, F. Wilhelm, R. Schouten, C. Har- 
mans, T. Orlando, S. Lloyd, and J. Mooij, Science 290, 773 

(2000) . 

[3] D. Vion, A. Aassime, A. Cottet, P. Joyez, H. Pothier, C. Urbina, 
D. Esteve, and M. H. Devoret, Science 296, 886 (2002). 

[4] Y S. Greenberg, A. Izmalkov, M. Grajcar, E. Il'ichev, 
W. Krech, and H.-G. Meyer, Phys. Rev. B 66, 224511 (2002). 

[5] Y S. Greenberg, A. Izmalkov, M. Grajcar, E. Il'ichev, 
W. Krech, H.-G. Meyer, M. H. S. Amin, and A. Maassen 
van den Brink, Phys. Rev. B 66, 214525 (2002). 

[6] E. Il'ichev, N. Oukhanski, A. Izmalkov, Th. Wagner, M. Gra- 
jcar, H.-G. Meyer, A.Y Smirnov, A. Maassen van den Brink, 
M. H. S. Amin, and A.M. Zagoskin, Phys. Rev. Lett. 91, 097906 
(2003). 

[7] D. Born, V. 1. Shnyrkov, W. Krech, T. Wagner, E. Il'ichev, 
M. Grajcar, U. Hubner, and H.-G. Meyer, Phys. Rev. B 70, 
1 80501 (R) (2004). 

[8] J. M. Raimond, M. Brune, and S. Haroche, Rev. Mod. Phys. 73, 
565 (2001). 

[9] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. J. 
Schoelkopf, Phys. Rev. A 69, 062320 (2004). 
[10] A. Wallraff, D. Schuster, A. Blais, L. Frunzio, R.-S. Huang, 
J. Majer, S. Kumar, S. M. Girvin, and R. J. Schoelkopf, Nature 
431, 162 (2004). 

[11] A. A. Abdumalikov, Jr., O. Astafiev, Y Nakamura, Y A. 

Pashkin, and J.S. Tsai, Phys. Rev. B 78, 180502(R) (2008). 
[12] R. Bianchetti, S. Filipp, M. Baur, J. M. Fink, M. Goppl, P J. 

Leek, L. Steffen, A. Blais, and A. Wallraff (2009), 0907.2549. 
[13] M. Grajcar, A. Izmalkov, E. Il'ichev, T. Wagner, N. Oukhanski, 

U. Hubner, T. May, 1. Zhilyaev, H. E. Hoenig, Y S. Greenberg, 

et al., Phys. Rev. B 69, 060501(R) (2004). 
[14] F. Yoshihara, K. Harrabi, A. O. Niskanen, Y Nakamura, and 

J. S. Tsai, Phys. Rev. Lett. 97, 167001 (2006). 
[15] J. Hauss, A. Fedorov, C. Hutter, A. Shnirman, and G. Schon, 

Phys. Rev. Lett. 100, 037003 (2008). 
[16] D. 1. Schuster, A. Wallraff, A. Blais, L. Frunzio, R.-S. Huang, 

J. Majer, S. M. Girvin, and R. J. Schoelkopf, Phys. Rev. Lett. 

94, 123602 (2005). 



